Vortices in attractive Bose-Einstein condensates in two dimensions 
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The form and stability of quantum vortices in Bose-Einstein condensates with attractive atomic 
interactions is elucidated. They appear as ring bright solitons, and are a generalization of the 
Townes soliton to nonzero winding number m. An infinite sequence of radially excited stationary 
states appear for each value of m, which are characterized by concentric matter-wave rings separated 
by nodes, in contrast to repulsive condensates, where no such set of states exists. It is shown 
that robustly stable as well as unstable regimes may be achieved in confined geometries, thereby 
suggesting that vortices and their radial excited states can be observed in experiments on attractive 
condensates in two dimensions. 
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The study of vortices has a long and illustrious scien- 
tific history reaching back to Helmholtz and Lord Kelvin 
in the nineteenth century Vortices associated with 
quantized circulation are a central feature of superflu- 
idity . Singly- 0, 01 and multiply-quantized [3 vor- 
tices have been observed in Bose-Einstein condensates 
(BECs) with repulsive atomic interactions. Complex vor- 
tex structures have been shown to be stable in repulsive 
BEC's, including vortex dipoles and vortex rings 0- 
The nonlinear Schrodinger equation (NLSE), which pro- 
vides an excellent description of BECs at the mean-field 
level , supports vortex solutions, which have been stud- 
ied extensively in the case of repulsive atomic interac- 
tions The main goal of this Letter is to clarify the 
meaning and nature of single vortices and their excited 
states in BECs with attractive interactions, and thus en- 
courage experimental exploration of stable and unstable 
two dimensional (2D) BECs. 

Solutions to the NLSE with attractive, or focusing non- 
linearity, in contrast to repulsive, or dcfocusing nonlin- 
earity, are unstable in free space in three dimensions (3D) 
and are stable in one dimension (ID) The imposi- 

tion of an external potential in 3D can produce a long- 
lived mctastable regime ^3 ^3 ■ the metastablc and 
unstable regimes, growth and collapse cycles and im- 
plosion have been studied in 3D. In the stabl e reg ime, 
bright soliton propagation [l^ and interactions jl7l | in a 
waveguide have been investigated in ID [l^ . 

The critical dimensionality for the NLSE is 2D [lllligt. 
We will show that quantum vortices and their radially 
excited states can be made robustly stable in confined, 
attractive 2D condensates and are a generalization of 
the Townes soliton \2(\ to nonzero winding number m. 
The Townes soliton is fundamental to understandin g th e 
self-similar collapse of solutions to the 2D NLSE [21|. 
Its generalization to winding number |m| > 1 has been 
studied in the context of optics, where such solutions are 
called "ri ng-p rofile solitary waves" or "spinning bright 
solitons" [231 . This is in fact the attractive analog of 



the well-known single vortex solution in repulsive con- 
densates 0], as we will show. An example contrasting 
the form of vortices in attractive and repulsive BECs is 
illustrated in Fig. ^a) and ^b) . 

In previous studies of quantum vortices in condensed 
matter and optical systems with attractive non- 
linearity, the phase variation derived wholly from circula- 
tion of matter about the central vortex core. In this work, 
we investigate the most general type of single- vortex sta- 
tionary solutions in attractive BECs for which the phase 
of the order parameter also alternates sign along radial 
lines. For each winding number m, we find a denumer- 
ably infinite set of radially excited states characterized 
by the successive formation of nodes at r = c». An ex- 
ample of such an excited state is illustrated in Fig.^c). 
In contrast. Fig. ^d) shows how extremely different is 
the repulsive case. 

Consider an order parameter of the form 

V'(r, 0, t) ^ f„i{r) exp(im(^) exp{-ifit/h) exp(i6'o) , (1) 

which solves the 2D NLSE 

[-(nV2M)V2 ^ -I- gaolV'P + V{r)]ij = zfidti^ , (2) 

with V{r) a central potential in two dimensions, m 
the winding number, M the atomic mass, 320 = 
j ^[mZj^J^t^ the 2D atomic interaction 
strength, < the s-wave scattering length, and /i 
the chemical potential or eigenvalue. In Eq. (O it was 
assumed that the BEC remains in the ground state in 
a harmonic trap of angular frequency in the z direc- 
tion |23|, so that r = ^ +2/^. Then, taking 0o = and 
defining 77„(x) = \g-m I lA! 'm{r) , x = {\/'^M\^i\/h)r, 
one obtains a rescalcd 2D NLSE of form 
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+ ^?,n-"^(x)^7m+W>-0' (3) 



where ~ sgn(/i) = ±1. The solutions to this nonlin- 
ear second order ordinary differential equation describe 
quantum vortices and their radially excited states in an 
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FIG. 1; A quantum vortex of winding number m = 1: (a) 
attractive case; (b) repulsive case. A radially excited state: 
(c) the first excited state in the attractive case; (d) in the re- 
pulsive case, a radially excited state requires an infinite num- 
ber of nodes and asymptotically resembles the Coulomb func- 
tion |2^ . The radial dependence of the order parameter of an 
infinitely extended condensate is depicted. Note that all axes 
are dimensionless. 

attractive BEC in an external potential V{x). Examples 
for = —1 and V ^ are shown in Fig.^a) and^c). 

The form of the radial wavcfunction rjm can be ob- 
tained from Eq. ^ by nmnerical shooting methods |25l |. 
This requires the initial conditions rjmixo), v'mixo): 
which can be obtained via a power series around x = 0: 

where the aj are coefficients. Note that i]m{x) ^ aoX™ 
as X ^ 0. Upon substitution into Eq. ||3Jl and solution of 
the resulting simultaneous equations, one finds that all 
coefficients aj for j 7^ can be expressed as a polynomial 
in powers of gq. Thus the power series shows that ag, the 
coefficient of the first nonzero term in the series, together 
with the winding number m, is sufficient to determine the 
entire solution. Note that we take ao > and m > 0; 
solutions identical in form exist for oq < and/or to < 0. 

We first consider the case of no external potential, 
^(x) = 0- By following the form of the wavcfunction 
as a function of oq, one observes its entire development. 
The case of to = 1 and (7^ = — 1 is illustrated in Fig. |21 
For very small og, the linear, divergent Bessel function 
solution Krn{x) is recovered for small rjmix)- However, 
the nonlinearity regulates the divergence as 7y,„ becomes 
on the order of unity. It subsequently oscillates around 
unity, with the oscillations damping away as x ^ 00 
(Fig. I3a)). As ao approaches a critical value a™''*°'^, the 
oscillations are pushed out towards x = c» and a local- 
ized central peak appears in the wavefunction near the 
origin (Fig.[3^b)). For uq ~ a™*'''^ precisely a node 
forms at X = 00 and a quantum vortex is obtained. As 
flo is increased, the node moves inwards and oscillations 
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FIG. 2: The development of excited states of a quantum vor- 
tex in an attractive BEC. (a) For small ao, no quantum vortex 
is formed, (b) As ag a™'^''^'' from below (solid black curve) 
or above (dashed red curve), the oscillations are pushed out 
towards x ~ oo- For oq = a™'^''"' a node appears at x = 00, 
and a true vortex appears, (c) For still larger ao, the node 
moves inwards, (d) As ao approaches a second critical value, 
the oscillations are again pushed out, until a second node ap- 
pears at X = 00. In this way one obtains the first excited state. 
An infinite sequence of excited states can be constructed in 
this way. Note that the case of an infinitely extended conden- 
sate of winding number m = 1 is depicted and all axes are 
dimensionless. 



resume beyond it (Fig.l^Ja)). To move the oscillations to- 
wards infinity again, one must increase og towards a sec- 
ond critical value (Fig.|2Id)). When this value is reached 
precisely, one obtains the first excited state, with no os- 
cillations past the first node, and a second node appears 
at X = oo- In this way one can construct an infinite set 
of excited states by increasing Oq to successive critical 
points. These critical points are always characterized by 
the formation of an additional node at x = 00. 

The appearance of a denumerably infinite set of crit- 
ical points is in strong contrast to the case of a vortex 
in a repulsive condensate. As we have shown in previous 
work [221, for a repulsive BEC there is only one critical 
value of flo, i.e., 0™''*°''. Larger values of oq result in an 
infinite number of nodes and the wavefunction asymptot- 
ically resembles the Coulomb function, as illustrated in 
Fig. ^d). We called these ring solutions, in contrast to 
vortex solutions, as they have a different analytic struc- 
ture and asymptotic behavior. Thus, in an infinitely ex- 
tended system, vortices in repulsive condensates cannot 
be radially excited in a stationary way, whereas in attrac- 
tive condensates they have a denumerably infinite set of 
excited states. We note that, for positive chemical po- 
tential, i.e., CT^ = +1, one can also find ring solutions in 
attractive condensates. 

For winding number m = and oq = 2.20620086 . . ., 
one obtains the Townes soliton. Increasing ag results in 
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the formation of successive nodes at x = oo, just as for 
m = 1, i.e., one finds the radiaUy excited states of the 
Townes sohton. In this special case, it has been formally 
proven that an infinite set of radially excited states cor- 
responding to the formation of nodes exists, and that the 
Townes soliton is the unique "ground state" in that se- 
quence . It is in this sense that the vortex solutions we 
have described are generalizations of the Townes soliton 
to nonzero winding number. All such solutions wherein 
a node has formed at infinity are normalizable. 

In order to study attractive BECs in experiments, it is 
vital to consider stability properties in confined systems. 
The special stability properties of 2D can be illustrated 
by a simple variational study. Consider the variational 
ansatz 



(5) 



where tq and A are variational parameters, subject to 
the power law potential V(r) = Vqt^ , j > 0. Integrating 
Eq. over 'ip*{r, 4>, t), one finds simultaneous equations 
for the chemical potential and (720- Then, using the nor- 
malization j d^r|-!/)p = TV to eliminate where N is the 
total number of atoms, one obtains parametrically 
in ro: 
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where J\f = MN\g2i)\/2TTn = 2\as\N Mluz/ V^nh is the 
dimensionless nonlinearity. The solution is radially sta- 
ble when the Vakhitov-Kolokolov (VK) criterion [2^ 
d^/dM < is satisfied. One finds that the VK crite- 
rion always holds for M < Mc-, where 



Nc = 2V^r(m + 2)/r(TO -fl/2) . 



(8) 



Equation (|SJ) is independent of both Vq and j, i.e., ra- 
dial stability docs not depend on the details of any 
positive power law potential. In the limit in which 
i — > 00, one obtains the same result for a cylindrical 
hard- wall potential. When M > A/'c, the total energy 
E[^] = N[n - (52d/2) / d'riVH -> -00 as ro ^ 0, 
meaning that the wavefunction is unstable and implodes. 
When M < A/'c, the energy has a global minimum at some 
finite tq. In contrast, in 3D, the minimum, when it exists, 
is alwa ys l ocal, so that, at best, one obtains metastabil- 
ity mil. 

We now consider in detail the most experimentally rel- 
evant case, i.e., j = 2 and Vb = Muj"^ a harmonic 
trap. Note that typical trapping frequencies are on the 
order of = 27r x 100 Hertz, which gives a time scale of 
T = 2'k/uj = 10 ms. For no external potential, V{r) ~ 0, 
it is well known that the Townes soliton (m = 0) is radi- 
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FIG. 3: (color online) Stable regimes of quantum vortices in 
confined attractive BECs. Shown are the chemical potential 
spectra and linear instability frequencies as a function of the 
nonlinearity for (a)-(b) the first excited state of the Townes 
soliton, and (c)-(d) the first excited state of a singly quantized 
vortex m = 1, all in a harmonic potential. In (b),(d) the 
winding number of the instability mode is denoted by q and 
the instability time given by 2n /lm{Qq). Note that all axes 
are scaled to harmonic oscillator units. 



ally unstable, while the quantum vortex of winding num- 
ber TO = 1 (Fig.^a)) is azimuthally unstable. The varia- 
tional analysis of Eqs. 10-® suggests that the addition 
of a confining potential stabilizes the solutions radially 
for general m. However, this simple variational study 
did not allow for azimuthal stability. Both radial and 
azimuthal stability can be determined by linear stabil- 
ity analysis, i.e., the Boguliubov-de Gcnnes equations 
(BDGE) yjj], a standard method which we do not re- 
produce here, for the sake of brevity. The case of gen- 
eral m without radial excitations has been studied previ- 
ously [iil Isof . It was found that for sufficiently small 
JV the vortex with to = 1 is stable, while for m > 2 the 
solution is unstable to quadrupole oscillations, though 
instability times can be much greater than T. We can 
therefore say that for small f\f vortices of winding num- 
ber TO > 2 arc experimentally stable. 

Applying the BDGE's, we studied the azimuthal and 
radial stability of the radially excited states of the Townes 
soliton. The results for the first excited state, which has 
a matter-wave ring separated from a soliton core by a 
radial node, are illustrated in Fig. El;a)-(b). The wind- 
ing number of each Boguliubov mode is denoted by q. 
The frequency of each mode is denoted by flq with in- 
stability time tq = 2ir /lm{ilq) . From Fig. OJb), it is 
apparent that a radial instability (q ~ 0) occurs for all 
A/", since Im(r2o) ^ 0. However, one finds experimen- 
tal stability for small M, since tq/T ^ 1. Other modes 
become unstable at higher JV, starting with the dipole, 
q = 1. For large Af the dominant instability occurs in 
the quadrupole mode. 
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In Fig. |2Ic)-(d) is shown the same stabihty analysis 
for the first excited state of the m = 1 vortex, which has 
two concentric matter-wave rings separated by a radial 
node. Unlike the excited states of the Townes soliton, 
here the solution is linearly stable for <5.5. The solu- 
tions continues to be experimentally stable up to ~10, 
although a dipole instabihty occurs for ti » T. For 
large Af the dominant instability is radial. For the ra- 
dially excited states of both the Townes soliton and the 
m = 1 vortex, the VK criterion fails, as is apparent in 
Figs. |2Ja) and|3{c). This is first instance we know of 
where this criterion, which so far as we know has never 
been proven formally, fails. For winding number m ~ 2 
we find that the first excited state is linearly most unsta- 
ble to quadrupolar excitations, although for A/" < 10 it is 
experimentally stable. 

In experiments, one may access the stable regime of 
quantum vortices and their excited states via a Feshbach 
resonance, which allows for extremely small scattering 
lengths on the order of the Bohr radius; this technique 
was used to successfully create bright solitons 0, Il7| . 
Vortices may be created by rotating the condensate; al- 
ternatively, one may start with a single vortex in a repul- 
sive condensate and adiabatically switch the scattering 
length from positive to negative. Excited states may be 
made by phase imprinting through a pinhole mask, by 
passing a tightly focused laser pulse through the conden- 
sate center, or by utilizing a doughnut mode of a laser. 
It seems likely that the most useful approach will use 
these methods to first create dark ring solitons j^^l on a 
repulsive condensate with a central vortex |2.'^ , and then 
tuning the scattering length to be negative. An impor- 
tant point in the assumption of the 2D regime is that, 
in order for our stability criterion to be valid, the con- 
densate must remain oblate, so that excitations in z do 
not occur; this may be achieved by a correct choice of 
experimental parameters. 

In conclusion, we have shown that quantum vortices 
and their radially excited states in attractive BECs can 
be created stably in confined systems. We contrasted 
vortices in attractive BECs, which can be thought of as 
ring bright solitons or spinning Townes solitons, to their 
counterparts in repulsive BECs. We showed that there 
exists a denumerably infinite set of excited states which, 
in an infinitely extended condensate, correspond to the 
creation of nodes at r = oo. In a harmonic trap, these 
are stable or experimentally stable for sufficiently small 
nonlinearity. In contrast to 3D, there is no metastability. 
We note that, in an optics context, vortices, or "spinning 
bright solitons," can also be stablized by competing non- 
linearities, i.e., a defocusin g qu intic nonlinearity with a 
focusing cubic nonlinearity [22|- 
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